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The ground-state energy of three-body ions {M^,M'^ 



evolves when the like-charge con- 
7n~) with the average of 



stituents are given different masses. The comparison of {mf, 

{mf,mf,m~) and {m2 ,m2 ,m~) reveals a competition between the symmetric term and the anti- 
symmetric one. The former dominates in the Born-Oppenheimer regime such as the (p,t,e) case, 
while the latter wins for H~-like systems with two negative light particles surrounding a heavy 
nucleus. A comparison is also made with the case of baryons in simple quark models with flavour 
independence. 
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I. INTRODUCTION 



There are many examples of molecules built out of con- 
stituents with identical charge which might acquire differ- 
ent masses when isotopes become involved. For instance, 
the H2 molecule is seen in several variants such as DT 
where deuterium and tritium replace the protons. In the 
field of exotic atoms and molecules, many configurations 
can be envisaged, such as (p, /i~,7r~) where two of the 
constituents have slightly different masses. 

The aim of this paper is to analyze the variation of en- 
ergy in the simplest case of the ground state of H^-like 
configurations. There is already a considerable literature 
on such systems {jyj, and tables of very accurate varia- 
tional energies are available. The mass-dependence of the 
energy has also been analyzed, and useful approximate 
formulas have been proposed. Here, we shall restrict our- 
selves to the specific problem of symmetry breaking. 

For consistency, the mass-dependence of the energy 
E{mi,m2) of the the {mj ,m2 ,Tn~) ion will be analyzed 
from the published (or recalculated) values at several 
neighbouring points, or from an estimate of the first and 
the second-order perturbation terms calculated with the 
wave function of a central configuration (M+, M"*", m~). 

The behaviour of the {mi,m2,m) systems as a func- 
tion of TOi and 7712 for a given third mass m and a given 
interaction potential was debated some years ago in the 
framework of confining quark models with flavour inde- 
pendence, that is to say, an interaction that does not 
change when the mass of any of the quarks is varied, 
in the same way as the Coulomb potential governs H~, 
Hj", or Ps~. This note offers the opportunity to make a 
comparison between the atomic and hadronic cases. 
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FIG. 1: Ground-state energy of (A/^, M^ ,m ) as a function 
of m/M, for m = 1. 



II. GENERAL CONSIDERATIONS 

The non-relativistic Hamiltonian is written, with an 
obvious notation, as 



H = 



P? 



P2 



2toi 27T12 



-El 

2m 



V . 



(1) 



where the potential V is the familiar V 



12 



The energy of any level is an increasing function of 



each inverse mass 



This is due to p^ being positive. 
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For the ground state, the energy is a concave function 
of any parameter entering the Hamiltonian linearly j2], 
in particular an inverse mass. The result is true for the 
ground state of any given angular momentum, or for the 
sum of the n first levels. 

In Fig. n the energy of the ground state of the sym- 
metric case, corresponding to the Hamiltonian 

^s(M,m) = ^[p?+p^] + ^p2 + F, (2) 

is shown as a function of the ratio x — m/M , for r?T, = 1. 
The monotonic and concave behaviour is clearly seen. 



The values of Fig. ^ can be taken from Ref. |^, but we 
have preferred to recalculate the energies using a similar 
variational method employing randomly chosen complex 
exponents }4J. The regular pattern indicates that a safe 
interpolation of the binding energy can be done from a 
few values of the masses. 

The concavity property in Fig. ^ shows that, within 
symmetric configurations, the state {AI~^ , M'^ ,m^) has 
its energy above the average of {mf ,mf ^m") and 
{m2 ,m2 ,'iTi~) if its inverse mass is taken the exact av- 
erage, namely, 



2M"i = m.-^ 



(3) 



The system of interest is, however, (mj", TO2 , f7z~), and 
it receives a downward shift with respect to its sym- 
metrized version (M+, M+,r7i^). Indeed, the Hamilto- 
nian (Q can be rewritten as 
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FIG. 2: Binding energy, Eb, of (Af+, M"'",m~) as a function 
of m/M for m = 1. 



H{mi, ?Ti2, rn) — Hs{M, m) + - 



-{pI-pD- (4) 



The symmetry-breaking term oc {m^ — m2 ) lowers the 
ground-state energy, as can be seen by applying the vari- 
ational principle to 10} , with the symmetric ground state 
of Hg as trial wave function. 

In short, when comparing the average of {m^ , mj , m~ ) 
and {m2 ,1712 ,m^) to the mixed state {mj ,m2 ,m^), 
there is a conflict between the concave behaviour in the 
average inverse mass of the positive particles, and the 
downward shift due to symmetry breaking. 

Starting from the configuration (M+ , Af + , m~ ) , one 
can write the expansions 



E{mi,m2) 



E{M,M) + Kei,s + Xh2.^s + --- (5) 

E{M, M) + Aei., -I- Pes,. + /i^ea.a + • • ■ 



where the common third mass is omitted and 

A, = (2m,)-i-(2A/)-\ A = (Ai+A2)/2, 
^ = (4mi)-i-(4m2)"' , (6) 

which is reduced to A = 0, Ai = — A2 = ^, if M^^ is 
taken as the exact average of inverse masses mj~ and 
TO^ . In the latter case these expansions reduce to 



E{mi,m2) = E{M,M)+fi'^e2,a 
E{mi,mi) + E{m2, TO2) 



= EiM, M) + ^J,'e2,s + 



(7) 



The coefficients e can be estimated cither by fitting 
the energies computed for the neighbouring values of the 
constituent masses, or estimated from the wave function 
of {M^ ,M'^,m~), using perturbation theory. The gen- 
eral results listed above translate into ei . > for each 
level, and £2,. < and £2,0 < for the ground state in 
any sector of conserved quantum numbers. 



III. RESULTS 



A. Stability 



It can be shown analytically [3 that (A/+, A'/+, m^) is 
stable against dissociation into (M+, m~) and an isolated 
M"*", for any value of the mass ratio M/m. The binding 
energy 



Eb = E{M+,M^ 



-E,,,{M+,m-) 



(8) 



improves dramatically as M increases, in the region 
where M > m. However, Eb/Ei-h is not monotonic in 
the full range of values of M/m. In Fig. 2, it is seen that 
the binding energy reaches its minimum at Af « 0.5 and 
then starts to increase again when the system goes to the 
H" limit. 



B. Concavity vs. symmetry breaking 

In order to analyze the behaviour of E(mi,m2,m) 
around the central configuration E{M, M, m) (see 
Eq. |SJ|) we introduce the two operators which determine 
the coefficients of expansion in ((Sj) , namely. 



Vs=p\+pI. and Va=V>\- V>1 
The perturbation theory gives 

e2,a = (*o|K(^o-ifo)-^K|*o), 



(9) 



(10) 



where Q — I — |^o)(^o| is the projection operator into 
the subspace orthogonal to |\l/o), and Hq, the symmetric 
Hamiltonian Hs{M,m) with m = 1. The first order 
perturbation for Va vanishes. 
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E 
-0.04774837464 
-0.08729940294 
-0.12069565567 
-0.14932813533 
-0.17418583295 
-0.19599515999 
-0.21530357938 
-0.23253249352 
-0.24801209057 
-0.26200507023 
-0.35268735314 
-0.40026703884 
-0.46608749725 
-0.50180015339 
-0.55024630859 
-0.57644241286 
-0.58833372341 
-0.59509329970 
-0.59836899224 



0.008647 
0.028959 
0.055480 
0.085138 
0.116148 
0.147453 
0.178428 
0.208710 
0.238095 
0.266477 
0.496963 
0.658243 
0.961169 
1.206036 
1.851059 
2.944942 
4.637278 
7.900453 
13.16189 



£2,8 

-0.001567 
-0.009631 
-0.025623 
-0.048893 
-0.078225 
-0.112387 
-0.150301 
-0.191087 
-0.234046 
-0.278636 
-0.757118 
-1.243673 
-2.754451 
-5.029963 
-20.99919 
-117.3044 
-595.6470 
-3598.177 
-18659.60 



£2, a 

-0.009206 
-0.059668 
-0.163192 
-0.314523 
-0.502118 
-0.713757 
-0.938895 
-1.169321 
-1.399056 
-1.623973 
-3.356412 
-4.319080 
-5.508756 
-6.049197 
-6.800187 
-7.767019 
-9.343615 
-12.52000 
-17.73284 
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FIG. 3: The ratio of the second order terms (symmetry vs. 
asymmetry) . 



N 


£2,s 


N 


£2, a 


600 10.453482 
800 10.453482 


-8940.3037 
-8940.3038 


600 
800 


-15.043224 
-15.043224 



TABLE I; The first and second order perturbation for the 
operators Va and Va for the ground state of the system 
{M+,M+,m~} with m" = 1. 



Results of numerical calculations are presented in Ta- 
ble I. From these data one can get, for instance, an ap- 
proximate value for the HD ion ground state, using data 
for Af = 3000 and Eq. ©, 



E{mp,md) = -0.597901 a.u. , 
to be compared with the exact solution 
E^^+ = -0.597897968 . . . , 



(11) 



(12) 



showing that this rough estimate provides about 6 sig- 
nificant digits for the energy. 

From these data one may define approximate bounds 
for the stability region by solving the equation. 



E{M, M) + fih2,a = -l(jj + ^f' + ^ 



(13) 



with respect to /i. Here we assume that mi > 7712 
and m = 1. The variation of fj, corresponds to mov- 
ing along an horizontal line in the triangle of normal- 
ized inverse masses |a|- The value of /i obtained from 
Eq. I|13(l for M = 0.1 is /i w —0.55, that corresponds 
to nil — 0.112 and m2 ~ 0.0901, the ratio of these two 
masses is mi/m2 — 1.25. 

The ratio of the symmetric to asymmetric second-order 
terms (i.e., concavity vs. symmetry breaking) is shown in 
Fig. 3. In the adiabatic region M 3> m, the symmetric 
term dominates. An inverse hierarchy is observed in the 



TABLE IL The first and second order perturbation for the 
operators Vs and Va for the ground Isag state of the H^. 



atomic regime M <^ m, with the asymmetric term dom- 
inating by about one order of magnitude, implying an 
inequality 



2E{mi,m2) < E{mi,mi) + E{m2,m2) 



(14) 



C. H2 "gerade" and "ungerade" ground states 

To get a different insight to the problem, we now com- 
pare the concavity and symmetry breaking terms for the 
H2 molecular ion ground states of different symmetry, 
starting from the molecule with identical nuclei. The 
2pau (or "ungerade") state is bound by a very shallow 
potential and lies very close to the threshold. On con- 
trary, the IscTg potential well supports the existence of 
19 vibrational states. Hence it is anticipated that open- 
ing the coupling between the "gerade" and "ungerade" 
states lowers the ground 2p(7„ state more rapidly than 
the ground Isag state. 

Numerical results are presented in Tables II and III. 
The "ungerade" 62,0 exhibits much slower convergence, 
since the intermediate state, or the first-order perturba- 
tion wave-function, a symmetric function being solution 
for the Isag potential at about zero energy, has multiple 
nodes that the variational approximation should repro- 
duce with a reasonable accuracy. This makes the numer- 
ical calculations more delicate. 

The results in the Tables II and HI illustrate our con- 
jecture that that the symmetry breaking effect for the 
case of 2pcr,j state is much stronger than for the "ger- 



N 



£2., 



N 



£2,0 



600 1.0497471 
800 1.0497471 
1000 1.0497471 



-99.715284 
-99.725303 
-99.725735 



1000 -40051.369 
1400 -40342.387 
1800 -40376.809 



TABLE III: The first and second order perturbation for the 
operators Vs and Va for the ground 2pau state of the H J . 



It was even conjectured that Eq. (|16|l might be gener- 
ally true, though resisting proof. Then Lieb found 
counterexamples, which were further extended in [ 
Lieb also found wide classes of potentials for which the 
inequality (|16|l is rigorously proven. Note that the quark 
problem does not include the case of H^-like ions, since 
the interquark potential V = '^v{rij) is assumed to be 
fully symmetric. 



ade" ground state, while the concavity term is almost 
negligible. 



IV. COMPARISON WITH BARYONS 

In Quantum Chromodynamics (QCD), quarks are cou- 
pled to gluons through their colour. The static potential, 
which serves as first approximation to describe the inter- 
action, is thus independent of the mass of the quarks 
experiencing this potential. This property, called flavour 
independence, stimulated several studies on how the spec- 
trum evolves in a given confining potential, as a function 
of the quark masses. For references, see Refs. ^^■ 

In the case of mesons, once the centre-of-mass motion 
is removed, the interaction Hamiltonian depends only 
upon the reduced mass fi, and the inequality [3, Q 

2i?2(mi,m2) > E2{mi,mi) + £2(1712,1712) , (15) 

holds whatever flavour-independent potential V{ri2) is 
assumed. 

In the case of baryons, it was always found for plausible 
flavour- independent potentials [9j, that 

2E^{mi,m2,7n) > E^{mi,mi,'m)+E^{m2,m2,7n) . (16) 

The approximate equality observed for experimental 
baryon masses (the contribution of constituent masses 
is identical on both sides) is due to a cancelation of this 
concave behaviour with the convex behaviour of hyper- 
fine corrections. 



V. CONCLUSION 

In this article, we have studied the effect of break- 
ing the TTT-i = 77i2 symmetry in the ground state 
(m^,mj,?7i~). Two regimes are identified. For mi and 
7712 , much heavier than the negatively charged particle of 
mass 777, the heavy particles experience an effective two- 
body potential. Then the dynamics depends mainly on 
the sum of inverse masses, a property that becomes exact 
in the case of two-body systems, or in the case of mesons 
in a flavour-independent quark-antiquark potential. In 
this situation, the inequality 116(1 is observed. Now, for 
77ii and 7772 of the order of 777 or smaller, in the H^ limit, 
the asymmetric term (Pi— P2) plays a more important 
role, and the reversed inequality is observed. 

For the adiabatic systems, which have both symmetric 
and antisymmetric bound states, the symmetry breaking 
effect is much stronger for the upper antisymmetric state, 
and is the dominant contribution, while the concavity 
term (p^-|-p|) is rather negligible. 

We intend to study the case of excited states, and more 
complicated structures, such as hydrogen-moleculc-like 
,m^,m~), and to extend the results of 



states (777^, 

our present work to the higher order perturbation terms. 
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